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ABSTRACT 

A list of 205 7-ray strong objects was reported recently as a result of a 3-month integration with the Large 
Area Telescope on board the Fermi Gamma-Ray Space Telescope. We attempted identification of these objects, 
cross-correlating the 7-ray positions with VLBI positions of a large all-sky sample of extragalactic radio sources 
selected on the basis of their parsec-scale flux density. The original associations reported by the Fermi team 
are confirmed and six new identifications are suggested. A Monte-Carlo analysis shows that the fraction of 
chance associations in our analysis is less than 5 per cent, and confirms that the vast majority of 7-ray bright 
extragalactic sources are radio loud blazars with strong parsec-scale jets. A correlation between the parsec- 
scale radio and 7-ray flux is supported by our analysis of a complete VLBI flux-density-limited sample of 
extragalactic jets. The effectiveness of using a VLBI catalog to find associations between 7-ray detections and 
compact extragalactic radio sources, especially near the Galactic plane, is demonstrated. It is suggested that 
VLBI catalogs should be used for future identification of Fermi LAT objects. 
Subject headings: radio continuum: galaxies — galaxies: active — catalogs 



1. INTRODUCTION 

The Fermi Gamma-Ray Space Telescope (previously 
known as GLAST) was successfully launched in June 2008. 
Even with the mode rn Large Area Telescope on board (LAT, 
Atwood et al. 2009), the positional uncertainty of 7-ray mea- 
surements in the energy range from about 20 MeV to more 
than 300 GeV remains relatively poor, typically 3' to 20'. 
This provides a challenge for 7-ray source identification. A 
large fraction of the high energy 7-ray sources detected by the 
EGRE T telescope on board the Compton Gamma Ray Obser- 
vatory dHartman et aljfl999l) was identified with blazar s (e.g., 
iMattox et alj200UISowards-Emmerd et alj2003Ll2004l) . This 
fact critically helps in identifying objects detected by the 
Fermi LAT. In order to help the identification process, several 
sample s of blazars and blazar candidates were constru cted re- 
cently dHealey et al]l2007l 120081: iMassaro et al.ll2009h . These 
compilations of many thousands of extragalactic objects, to- 
gether with other catalogs covering radio to 7-ray bands, were 
successfully used fo r identifying the LAT-detected sources by 
lAbdo et al.1 (l2009bh . 

Results of parsec-scale Very Long Baseline I nterferome- 
try (VL BI) measurements were not applied by lAbdo et ail 
(2009b) in the process of identifying bright LAT sources. 
However, this can provide important extra information and 
improve the estimation of the probability of correct identifi- 
cation since VLBI filters out objects which do not host strong 
compact jets at parsec scales. The latter is a strong character- 
istic of a radio loud blazar. VLBI can especially help in the re- 
gion around the Galactic plane where the available multi-band 
coverage of the sky is much poorer due to Galactic absorption. 
An obvious weakness of the VLBI approach is the fact that 
radio weak objects (e.g., high energy selected BL Lacs) are 
missed in the currently available large VLBI all sky surveys. 

During the last several decades, a number of large VLBI 
surveys were conducted covering the frequency range be- 
tween 2 and 100 GHz (Beas lev et all [2002; Fomalo nt et alj 



2000; iPollack et al] l200l iKovalev et al J 120051 lOiha et al 



2004 



Lovell et alj|2004t iHelmboldt et all 120071 lLanvi et al 



20051: Lee et al.l 120081 e. g.,) including ded icated surveys in 
the Galactic plane (e.g., iPetrov etlil]l200 7f). The larj 



rgest in 

terms of the number of sources observed is a VLBI effort at 
2 and 8 GHz with a prime g oal to construct an extragalac- 
tic celestial reference frame dMa et al.lfl998l iFev et alj |2004) 
and to find extragalactic sources which are com pact at par- 
sec sc ales and are suitable for phase referencing (Petr ov et alj 
120081 and references therein). We note that while the m ajor- 
ity of those are of a blazar type, Helmboldt et al. (2007) have 
found 10 per cent of compact extragalactic radio sources to 
be optically associated with galaxies from the Sloan Digital 
Sky Survey. In this paper we use the 2/8 GHz VLBI surveys 
to associate the early Fermi LAT bright source list with ra- 
dio counterparts at parsec scales, make a statistical check on 
the reliability of these identifications, and analyze the radio 
properties of the VLBI counterparts of the 7-ray bright AGN. 

2. VLBI CATALOG IN USE 

The VLBI positions for more than 4,000 extragalactic 
sources covering the whole sky were used in this analy- 
sis 1 . They are based on 2 and 8 GHz VLBI observa- 
tions, most of which were performed in the period 1994 to 
2007 dMa et al 1 119981: IFev et alj 120041 iBeasjev etail 120021: 
i Fomalont et al.ll2003UPetrov et alJl2005ll2006t IKovalev et all 
l2007t IPetrov et al] 120081. 12009allbl) . The accuracy of VLBI 
astrometric position determinations is typically better than 
1 milliarcsecond. For details of the observati ons, data pro- 
cessin g and analysis see the cited papers. IKovalev et al] 
(2007) have made an additional dedicated effort to construct 
a complete parsec-scale flux-density-limited sample with the 
following characteristics: flat radio spectrum with spectral in- 
dex a > -0.5 (S ~ v a ), parsec-scale 8 GHz flux density (in- 
tegrated flux density over VLBI map) 5Vlbi > 0-2 Jy, decli- 
nation 5 > -30°, Galactic latitude \b\ > l.°5. Known Galactic 
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FIG. 1, — Left panel: Absolute difference AP bet ween Fermi LAT po sitions and astrometric VLBI positions for positively associated sources from Table [T] 
versus square root of Fermi LAT Test Statistics from Abdo et al. 1 2009b). Middle panel: Fermi-WhJil AP difference versus Fermi LAT 95 per cent confidence 
radius from Abdo et al. Right panel: Relative difference AP/r^ versus square root of Fermi LAT Test Statistics. 



objects were excluded. The full all-sky astrometric 2/8 GHz 
VLBI source (AVS) catalog of 4338 objects is used for cross- 
identification with the 3 month F ermi LAT bright source 
(LBS) list from lAbdo etalj (|2009b) while the complete flux- 
density-limited VLBI source (FVS) sample of 1848 objects is 
used for statistical tests and analysis. 

3. Fermi LAT - VLBI CROSS -IDENTIFICATION 

A number of authors have suggested a close connection 
between the bright 7-ray a nd parsec-scale radio emission 
from analyzing EGRET (e.g., Lahteenmaki & Valtaoia 2003) 
and early Fermi (e.g jKovalev et alJl2009tlLister et al.l l2009b) 
data. They concluded that the related 7-ray and parsec-scale 
radio emission should originate in spatially close regions. 
This means that VLBI positions might serve very well as an 
estimate of a location of the source of 7-ray emission detected 
by Fermi LAT for extragalactic objects. Taking this assump- 
tion as our starting point, we have applied a simp le approach 
to cros s-identify bright Fermi LAT sources from Abd o et all 
(2009b). A Fermi source will be called successfully iden- 
tified with a VLBI one if the difference between the Fermi 
LAT position and the VLBI position of the AVS sample is 
less than the 95 per cent confidence radius, rq$, estimated by 
lAbdo et al.lJ2009bh . 

3.1. Results of the cross-identification 

The results of this cross-identification can be found in Ta- 
ble [T] We have successfully identified 111 out of 205 Fermi 
LAT sources. We note th at the full list of the 205 bright 7- 
ray sources presented by lAbdo et all d2009bl) includes both 
Galactic and extragalactic objects. No single Fermi detection 
was identified with more than one VLBI source. Figures[T|[2] 
present the results of these identifications. The square root 
of the Fermi Test Statistics (VTS) provides an estimate of 
the signal-to-noise ratio (SNR) for every Fermi LAT detec- 
tion. As expected, the difference, AP, between the Fermi 
position and that of its VLBI identification becomes smaller 
with increasing \/TS (Figure [TJ. There are no identifications 
found with high Fermi LAT SNR and high relative Fermi- 
VLBI positional difference, AP/rgs- This could be explained 
by an over-esti mated systematic u ncertainty which was added 
in quadrature (I Abdo et al.ll2009bT) while calculating rg$. The 
systematic uncertainty dominates the r^ value for high-SNR 



Fermi sources becoming less dominant with decreasing Fermi 
TS. Abdo et al. reported this systematic uncertainty to be con- 
servative, they expect it to improve. Th e distribution of AP 
(Figure |2]i agrees well with the results of lAbdo et alj (|2009b, 
Figure 9), who discuss this distribution in detail. 

Cross-identification of the Fermi LBS versus AVS lists 
resulted in 111 identifications out of 4338 VLBI sources 
(2.6 per cent of the VLBI sample), while LBS versus FVS 
lists comparison gave 79 identifications out of 1848 VLBI 
sources (4.3 per cent). We note that the fraction of identifi- 
cations within the complete parsec-scale flux-density-limited 
FVS list is about two times higher than that for the deeper 
AVS list of sources. This could be the case if the bright 7-ray 
sources from the LBS list "prefer" stronger parsec-scale radio 
counterparts, i.e. if there is a correlation between the 7-ray 
flux and parsec-scale radio flux density. We note that the flux 
density cutoff of the FVS sample is significantly higher than 
that for the AVS sample (200 mJy versus about 50 mJy). The 
correlation suggested is indeed confirmed by an analysis in 

§SI 

Comparison wit h the association results presented by 
lAbdo et al.l d2009bl) has shown the following. 105 Fermi LAT 
detections have identifications common to the Abdo et al. 
analysis and the current study. All the extragalactic associa- 
tions suggested by Abdo et al. but not found by us, have radio 
flux densities less than 0.2 Jy or declinations less than —30°, 
i.e. belong to the flux density and/or declination range where 
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FIG. 2. — Distribution of AP from Figure [T] Left panel: Absolute differ- 
ence AP. Right panel: Relative difference AP/r^. 
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TABLE 1 

VLBI ASSOCIATIONS OF THE Fermi LAT BRIGHT SOURCES 



LBS Name 


VLBI name 


AP 


''95 


VLBI R.A. 


VLBI Dec. 


SVLBI 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



OFGL J0643.2+0858 0640+090 0.043 

0FGLJ1 123.0-6416 1121-640 0.039 

OFGL J1328. 8-5604 1325-558 0.058 

OFGL J1604.0-4904 1600-489 0.037 

OFGL J1830.3+0617 1827+062 0.067 

OFGL J2001. 0+4352 1959+437 0.025 



0.121 06:43:26.4450 +08:57:38.0126 0.61 

0.125 11:23:19.4155 -64:17:36.1970 

0.142 13:29:01.1449 -56:08:02.6657 

0.078 16:03:50.6846 -49:04:05.5050 

0.097 18:30:05.9399 +06:19:15.9521 0.41 

0.069 20:01:12.8738 +43:52:52.8391 0.07 



NOTE. — TablefTlis presented in its entirety in the electronic edition of the Astrophysical Journal Letters. Here 
we present six new suggested identifications for convenience and guidance regarding Table's form and contents. 
Columns are as follows: (1) LAT bright source list catalog name lAbdo et aU2009bl) ; (2) IAU name (B 1950.0) of 
VLBI association; (3) Differe nce between Ferm i and VLBI position (deg); (4) Fermi LAT radius of 95 per cent 
confidence region (deg) from i Abdo et al. 2009b) is shown for comparison with AP; (5) J2000.0 Right Ascension 
of VLBI association in hr, min, sec; (6) J2000.0 Declination of VLBI association in angular hr, min, sec; (7) Total 
parsec-scale 8 GHz flux density of VLBI association object (Jy). 



our VLBI catalog is not complete. All the objects which are 
identified by us but have no counterparts suggested by Abdo 
et al. are located close to the Galactic plane — a region of 
special difficulty for identification by the standard method of 
Abdo et al., as discussed above. 

3.2. New identifications 

T he followin g six n ew identifications, not reported 
by lAbdo et alJ ([2009b), were found by our analy- 
sis: OFGL J0643.2+0858 (galactic latitude b = 2.°29), 
OFGL Jl 123.0-6416 (b = -3P02), OFGL J1328.8- 
5604 (b = 6.°41), OFGL J1604.0-4904 (b = 2? 54), 
OFGL J 1830.3+06 17 (b = 7.°54), OFGL J2001.0+4352 
(b = 1? 12). All these are within 10° from the Galactic 
plane. Details on these sources can be found in Table Q] 
We note that identifications with the same objects for two 
of these si x Fermi detec t ions were independently proposed 
recently bvlBassani et all (120091) for OFGL J2001 .0+4352 and 
iMirabal & Ha lpernl d2009b for O FGL J1830.3+0617 on the 
basis of a multi-band analysis. Mira bal & Halpernl (120091) 
have also discussed possible flat-spectrum radio source asso- 
ciations for OFGL J0643.2+0858 and OFG L J 1328 8 -5604 
which agree with our identifications. Mir aball (120091) 
presented an X-ray point-source catalog from Swift XRT ob- 
servations of unidentified Fermi sources including discussion 
of a potential counterpart for OFGL J 1604.0-4904. 

3.3. Chance association analysis 

We have performed a Monte-Carlo simulation to measure 
the chance coincidence probability. We kept the Fermi LBS 
catalog unchanged while the VLBI catalog was scrambled; 
every source in the VLBI catalog was shifted in a random po- 
sition angle for a random angular distance on the sky between 
0° and 5°. After that, cross-identification of the catalogs was 
done as described above. This exercise was performed 1,000 
times. The analysis was done independently for the LBS ver- 
sus AVS and LBS versus FVS lists. The mean number of 
chance associations was fou nd to be 4.5 (4.1 per cent of 1 1 1 
true identifications, see j33.U and 2.9 (3.7 per cent of 79 iden- 
tifications) respectively. This shows with a very high level of 
significance that (i) almost all of the associations found are 
firm identifications, (if) the original assumption of the method 
that the 7-ray sky is dominated by blazars with strong parsec- 
scale jets is confirmed. 

The described above method of VLBI catalog scrambling 
was applied to represent the non-uniformity of the AVS list, 



where sky regions of a special interest (e.g., the Galactic 
plane) are populated more densely. The 5° -radius was cho- 
sen as an optimal size to represent the non-uniformity being 
in the same time significantly greater than the r^ values. We 
note that, for a uniform sky coverage in an original catalog, 
this method reproduces uniformity in the Monte-Carlo test. If 
we consider as a match Fermi and VLBI sources with posi- 
tional difference less than rgs 12, we lose about 30 per cent of 
identifications (see Figures[T]l2} while the chance coincidence 
probability drops only down to about 2 per cent. 

4. VLBI PROPERTIES OF THE 7-RAY 
IDENTIFICATIONS FOUND 

We have performed a correlation analysis betwe en the av- 
erage F ermi LAT 100 MeV - 1 GeV photon flux (lAbdo et alJ 
2009b) versus 8 GHz parsec-scale flux density measured by 
VLBI between 1994 and 2008 (Figure[3]l. The non-parametric 
Kendall t test confirms a positive correlation in Figure [3] at a 
confidence level greater than 99.9 per cent for the 100 MeV - 
1 GeV energy band (the two upper limits were ignored by 
the test which did not affect the conclusion) . This agrees 
with re sults of a similar analysis performed bv lKovalev et alJ 
(2009) for simultaneous L AT 7-ray — 15 GHz VLBA MO- 
JAVE dLister et alj |2009a) measurements. Possible system- 
atics in this dependence resulting from the different proper- 
ties of two populations of extragalactic 7-ray sources (low 
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Fig. 3.— Average Fermi LAT 100 MeV - 1 GeV photon flux lAbdo et all 
2009b) versus 8 GHz parsec-scale flux density non-simultaneously measured 
by VLBA. The arrows represent an upper limit on the 7-ray flux. 
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strongly supports the correlation between 7-ray and parsec- 
scale radio emission , confirming again the recent finding by 
lKovalevetal](l2009b . 



5. SUMMARY 

It has been shown that VLBI provides a very efficient tool 
to identify bright 7-ray detections which have poor positional 
accuracy. Application of this method is especially important 
for Fermi LAT identifications near the Galactic plane where 
VLBI observations are affected by absorption and extended 
Galactic emission to a much lesser degree than other radio 
observations. It is suggested that this method be incorporated 
into the process of identification of the Fermi LAT catalogs 
and for estimation of systematics in 7 -ray positions. The 
results of 7-ray source identification by lAbdo et all ([2009b) 
have been confirmed. Six new identifications are reported for 
a follow up analysis, all within 10° from the Galactic plane. 
It is estimated that more than 95 per cent of the VLBI associ- 
ations found are firm identifications. 

Direct cross-correlation of the Fermi LAT and VLBI cat- 
alogs has confirmed with a very high level of s ignificance 
t he early finding of EGRET (lHartman et al.l ll999) and Fermi 
( Abdo etaLl 12009a) that most of the bright 7-ray sources 
on the sky are compact blazars. It was found that Fermi 
LAT preferentially "selects" the brightest objects from a flux- 
density-limited sample of radio loud parsec-scale jets which 
supports the hypothesis of a direct correlation and physical 
connection between emission in the 7-ray and radio bands. 



FIG. 4. — Distribution of parsec-scale flux density measured by VLBA at 8 
GHz for the complete flux-density-limited sample (FVS) with Svlbi > 0.2 Jy 
and declination S > —30°. Top panel: Positive Fermi- VLBI associations. 
Bottom panel: Fermi- VLBI non-associations. 

en ergy peaked versus high energy peaked) were mentioned 
by I Abdo et al.l (l2009al) . We do no t go into deta i led an alysis 
of this in the present paper since Kov alev et al.l ( 120091) have 
addressed this issue in their study. 

An even stronger test is done by analyzing the full FVS 
sample of 1848 sources which is complete down to 0.2 Jy, sep- 
arating them into LAT detected and LAT non-detected (Fig- 
ure @). We assume that LAT non-detected objects have sta- 
tistically lower 7-ray photon flux and therefore should have 
lower radio flux. The Kolmogorov-Smirnov test shows with 
a confidence greater that 99.99 per cent that the distributions 
of the parsec-scale flux density for the LAT detected and non- 
detected VLBI sources are drawn from different parent dis- 
tributions. The median flux densities of these two distribu- 
tions differ by a factor of 2.5: 0.84 Jy versus 0.34 Jy. This 
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